home *** CD-ROM | disk | FTP | other *** search
/ Language/OS - Multiplatform Resource Library / LANGUAGE OS.iso / a_utils / expanded.lha / expanded / src / Scalar.C < prev    next >
C/C++ Source or Header  |  1992-03-19  |  4KB  |  165 lines
  1. //
  2. // Linear-Affine-Projective Geometry Package
  3. //
  4. // Scalar.C
  5. //
  6. // $Header$
  7. //
  8. // William J.R. Longabaugh 
  9. // University of Washington
  10. //
  11. // Implementation of the linear-affine-projective geometry
  12. // package described in William J.R. Longabaugh, "An Expanded
  13. // System for Coordinate-Free Geometric Programming", Master's 
  14. // thesis, University of Washington, 1992.
  15. //
  16. // Copyright (c) 1992, William J.R. Longabaugh
  17. //   Copying, use and development for non-commercial purposes permitted.
  18. //   All rights for commercial use reserved.
  19. //   This software is unsupported and without warranty; it is
  20. //   provided "as is".
  21. //
  22. // ***********************************************************************
  23.  
  24. #include "Lap.h"
  25.  
  26. // ***********************************************************************
  27. // ***********************************************************************
  28. //
  29. // AugScalar Class
  30. //
  31. // ***********************************************************************
  32. // ***********************************************************************
  33. //
  34. // Public member functions
  35. //
  36. // ***********************************************************************
  37. //
  38. // Default augmented scalar equals 0.0:
  39. //
  40.  
  41. AugScalar::AugScalar(void) : m(1, 2)
  42. {
  43.   is_inf = FALSE;
  44.   m[0][0] = 0.0;
  45.   m[0][1] = 1.0;
  46. } 
  47.  
  48. // ***********************************************************************
  49. //
  50. // Create an augmented scalar with a finite value:
  51. //
  52.  
  53. AugScalar::AugScalar(Scalar v) : m(1, 2)
  54. {
  55.   is_inf = FALSE;
  56.   m[0][0] = v;
  57.   m[0][1] = 1.0;
  58. } 
  59.  
  60. // ***********************************************************************
  61. //
  62. // Create an augmented scalar with the value INFINITY:
  63. //
  64.  
  65. AugScalar::AugScalar(Infnum n) : m(1, 2)
  66. {
  67.   if (n == INFINITY) {
  68.     is_inf = TRUE;
  69.   } else {
  70.     errh.ErrorExit("AugScalar::AugScalar(Infnum)",
  71.            "Illegal value for building augmented scalar",
  72.                    ErrVal("Value of argument = ", n));
  73.   }
  74.   m[0][0] = 1.0;
  75.   m[0][1] = 0.0;
  76. }
  77.  
  78. // ***********************************************************************
  79. //
  80. // Automatically cast an augmented scalar into a scalar:
  81. //
  82.  
  83. AugScalar::operator Scalar()
  84. {
  85.   if (!is_inf) {
  86.     return (m[0][0] / m[0][1]);
  87.   } else {
  88.     errh.ErrorExit("AugScalar::operator Scalar()",
  89.            "Attempted to cast INFINITY to a scalar");
  90.   }
  91. }
  92.  
  93. // ***********************************************************************
  94. //
  95. // Get the value of the augmented scalar (This is just an explicit form of
  96. // casting to a scalar).
  97. //
  98.  
  99. Scalar AugScalar::Value(void)
  100. {
  101.   if (!is_inf) {
  102.     return (m[0][0] / m[0][1]);
  103.   } else {
  104.     errh.ErrorExit("Scalar AugScalar::Value(void)",
  105.            "INFINITE augmented scalar cannot return a value");
  106.   }
  107. }
  108.  
  109. // ***********************************************************************
  110. //
  111. // Need to do memberwise initialization, since Matrix class members need
  112. // to do memberwise initialization:
  113. //
  114.  
  115. AugScalar::AugScalar(AugScalar& a) : m(a.m)
  116. {
  117.   is_inf = a.is_inf;
  118. }
  119.  
  120. // ***********************************************************************
  121. //
  122. // Need to do memberwise assignment, since Matrix class members need
  123. // to do memberwise assignment:
  124. //
  125.  
  126. AugScalar& AugScalar::operator=(AugScalar& a)
  127. {
  128.   is_inf = a.is_inf;
  129.   m = a.m;
  130.   return (*this);
  131. }
  132.  
  133. // ***********************************************************************
  134.  
  135. void AugScalar::debug_out(ostream& c, int indent)
  136. {
  137.   char *ibloc = new char[indent + 1];
  138.   for (int i = 0; i < indent; i++) {
  139.     *(ibloc + i) = ' ';
  140.   }
  141.   *(ibloc + indent) = '\0';
  142.  
  143.   c << ibloc << ast;
  144.   c << ibloc << "Augmented scalar object\n";
  145.   if (is_inf) {
  146.     c << ibloc << "Value = INFINITY\n";
  147.   } else {
  148.     c << ibloc << "Value = " << (m[0][0] / m[0][1]) << "\n";
  149.   }
  150.   c << ibloc << ast;
  151.  
  152.   delete ibloc;
  153.   return;
  154. }
  155.  
  156. // ***********************************************************************
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.